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by 
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J.  P.  Godfrey 

ABSTRACT 

Goal  programming  has  become  a  popular  and  important  planning, 
tool  in  many  areas.  However,  the  main  direction  of  goal  programming 
research  has  been  in  formulating  models  instead  of  seeking  procedures 
that  would  provide  improved  computational  efficiency.  A  recent  paper  by 
Charnes,  Cooper,  Klingman,  and  Niehaus  has  exploited  the  specialized 
structures  of  a  class  of  convex  goal  programming  problems  to  obtain 
increases  in  efficiency.  In  this  paper,  an  analytical  formulation  for 
a  piecewise  linear  goal -interval  programming  problem  is  presented  and 
a  special  purpose  algorithm  is  developed.  The  algorithm  may  combine 
several  standard  simplex  pivots  into  a  single  pivot  through  the  use  of 
an  extended  minimum  ratio  test. 
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1.  Introduction 

For  a  substantial  period  of  time  following  the  development  of 
the  initial  goal  programming  model  for  a  scheme  of  executive  compensation 
by  Chames,  Cooper,  and  Ferguson  [11],  the  main  emphasis  of  goal  programming 
was  directed  towards  exploiting  this  tool  in  a  wide  variety  of  applications 
and  developing  formal  statistical  theory.  In  the  last  few  years,  there 
have  been  considerable  advances  in  the  "state  of  the  art."  With  the 
development  of  the  field  of  public  management  science,  more  sophisticated 
models  have  been  needed  to  solve  complex  problems  of  manpower  and  resource 
planning.  One  of  the  most  important  new  models  joins  goal  programming  and 
Markov  processes.  This  model  and  many  other  current  advances  can  be  found 
in  the  Office  of  Civilian  Manpower  Management,  U.  S.  Navy,  (OCMM)  series 
[14]. 

In  spite  of  the  wide  use  of  goal  prograimti ng ,  only  a  few  results 
have  been  obtained  that  will  secure  goal  programming  algorithms  that  are 
computationally  efficient.  A  paper  by  Chames,  Cooper,  Klingman,  and 
Niehaus  [12]  provides  a  specialized  solution  technique  by  a  transformation 
to  a  specialized  linear  interval  programming  problem  for  a  class  of  "goal 
functionals."  It  is  also  widely  known  that  the  problem  of  minimizing  the  sum 
of  the  absolute  values  of  the  goal  deviations,  the  L1  Norm  problem,  can  be 
placed  in  a  more  tractable  format  by  reduction  to  a  linear  programming 
equivalent.  In  the  case  of  solving  for  the  Least  Absolute  Value  Estimators 
(LAVE's),  the  special  structure  of  the  problem  has  enabled  the  creation  of 
special  purpose  primal  linear  programming  codes  [1,  2,  4,  5,  16  ]. 

The  difficulties  presented  by  nonlinear  convex  goal  functionals 
can  be  resolved  when  the  functional  is  separable.  In  this  case,  the  functional 
can  be  approximated  by  a  convex  piecewise  linear  functional  to  any  degree  of 


accuracy  that  may  be  desired  by  using  the  techniques  of  Charnes  and  Lemke 
[14].  The  piecewise  linear  functional  approach  was  used  by  Charnes  and 
Collomb  [8]  for  models  of  optimal  economic  stabilization  policy,  but  in¬ 
stead  of  using  the  classical  goal  programming  formulation  they  utilized 
the  concept  of  a  goal-interval  for  their  models.  Thus,  a  goal-interval 
functional,  in  a  single  variable  x,  will  assume  the  same  minimum  value 
over  an  interval  [g",  g+]  as  opposed  to  a  goal  programming  problem  which 
assumes  its  minimum  value  at  a  unique  point.  Collomb  [15]  also  studied 
goal-interval  models  for  decentralization  in  management  and  intertemporal 
analysis.  Charnes  and  Cooper  [9]  found  that  if  the  usual  goal  programming 
formulations  were  applied  to  assist  in  resource  allocation  decisions  for 
the  U.  S.  Coast  Guard's  Marine  Environmental  Protection  (MEP)  program, 
then  the  results  were  less  than  satisfactory.  Therefore,  to  obtain  more 
flexibility  in  the  models,  they  employed  the  goal-interval  approach  and 
obtained  results  more  in  line  with  the  desires  of  the  policy  makers. 

The  conceptual  flexibility  provided  by  goal-interval  programming 
would  be  negated  to  a  certain  degree  if  the  practical  implementation  for 
computational  purposes  caused  many  difficulties.  In  the  paper  by  Charnes, 
Cooper,  Klingman,  and  Niehaus  [12],  they  demonstrated  that  the  goal-interval 
programs  can  be  reduced  to  a  linear  programming  equivalent.  Hence,  the  full 
computational  and  interpretative  power  of  linear  programming  is  maintained. 
However,  since  there  are  many  variables  and  constraints  in  these  types  of 
problems,  it' would  be  very  advantageous  in  terms  of  money  and  computer  time 
to  have  an  efficient  algorithm  to  solve  the  problems.  A  specialized  primal 
simplex  algorithm  will  be  developed  to  exploit  the  convex  piecewise  linearity 
of  the  functional  for  deviations  not  in  the  goal  interval.  The  algorithm 


will  have  the  ability  to  combine  several  standard  simplex  pivots  into  a 
single  pivot  which  will  enhance  the  computational  aspect  of  applying  such 
techniques. 

2.  Analytical  Formulations  of  Goal-Interval  Programming 

A  desirable  aspect  of  goal-interval  model  is  that  it  presents 
a  mixture  of  a  satisficing  behavior  and  an  optimization  formulation.  The 
presence  of  separable  convex  piecewise  linear  functional  allows  the  con¬ 
venient  reduction  to  a  linear  programming  equivalent.  Collomb  [15]  pointed 
out  that  for  each  variable  there  are  thresholds,  which  define  regions  as 
'’satisfactory"  or  "less  than  satisfactory."  When  an  optimization  procedure 
is  applied  to  a  goal-interval  model,  then  the  optimal  solution  may  have  some 
variables  in  regions  which  are  "less  than  satisfactory."  Other  variables 
in  the  optimal  solution  may  be  permitted  to  vary  within  the  goal -interval 
since  this  provides  a  range  of  indifference  in  satisfaction. 

Figure  1  represents  a  goal -interval  model 
in  terms  of  a  single  incremental  variable,  and  it  is  similar  to  those  found 
in  Charnes  and  Collomb  [8]  and  Collomb  [15].  The  threshold  (goal)  levels 
are  given  by  the  g^'s,  and  the  goal-interval  is  [gj^,  g£].  The  w.  and  wt 
terms  are  the  relative  weights  (penalties)  assigned  to  the  "less  than 
satisfactory"  segments  of  the  functional. 

At  this  point  a  particular  type  of  problem  that  may  arise  in 
an  area  such  as  resource  allocation  is  formulated.  The  problem  for 
consideration  is  as  follows: 

m  m  +  + 

Minimize  t  w 76"  +  z  w.  6 
k=l  K  K  k=l  K  k 


(1) 


subject  to 


Sk  1  aklXl  *  ak2x2  * 
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+  + 
i  9,  + 
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b. 

i 


ihilxl  +  h12x2  * 


+  h.  x  <  b+ 
in  n  -  i 


i  =  1,  2, 


q 


+  - 

6k»  6k  1  0  for  k  =  1>  2,  . . . ,  m 


Rewritten  in  matrix  notation,  the  problem  becomes: 
MINIMIZE  w"  6‘  +  wV 

subject  to 

g~  -  6  <_  Ax  £  g+  +  6+ 

b  £  Hx  £  b+ 

6+,  6"  >  0 


The  system  b  <^Hx  <^b+  is  attached  to  represent  physical  constraints 
(limitations)  that  may  occur  and  that  must  be  satisfied  before  any  attempt  is 
made  to  resolve  the  goal  constraints.  The  goal-interval  for  the  k-th  goal 
constraint  is  [g^,  g*]  where  g^  <  g*.  For  the  case  of  gk  =  gk,  then  the  k-th 
goal  constraint  can  be  reformulated  as  a  usual  goal  equation  since  the  goal -interval 

a* 

has  been  reduced  to  a  unique  point.  If  a  vector  x  exists  such  that 
\x  >  ^k’  ^en  \  >  0  anc*  \x  =  9k  +  <5^ »  where  is  the  k-th  row  of  the 
matrix  A.  Thus  constraint  value  is  not  in  the  goal -interval ,  and  a  penalty  in 
the  amount  of  w£$*  is  incurred  for  the  excess  deviation  associated  with  this 
constraint.  A  simple  example  for  interpretation  of  this  goal -interval 
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constraint  may  arise  if  the  labor  force  is  divided  into  classes  where  a 
particular  class  may  be  associated  with  a  constraint  such  as  budgeting. 

If  the  operation  requires  more  employees  in  the  class  than  is  normally 
provided  for,  then  additional  employees  must  be  hired,  trained,  and,  in 
general,  "taken  care  of."  The  penalty  w£  associated  with  the  addition  of 
<$£  number  of  employees  may  arise  due  to  interest  on  loans  used  for  the 
expansion  or  due  to  whatever  other  criteria  are  used. 


The  previous  problem  (1)  is  one  of  the  simplest  formulations 
for  a  goal-interval  model.  The  assumption  of  constant  penalty  rates  for 


deviations  outside  of  the  goal -intervals  is  the  main  reason  for  this 
simplicity.  In  this  paper,  an  extension  of  problem  (1)  is  considered 
that  has  different  slopes  for  the  piecewise  linear  segments.  This  type 
of  problem  may  occur  is  when  the  class  of  employees  to  be  hired  are 
technicians,  and  an  additional  allocation  problem  for  the  existing  resources 
is  created.  If  a  large  number  of  technicians  must  be  hired  and  trained, 
then  there  are  limits  on  the  number  that  can  be  trained  efficiently  at  a 
time.  Thus,  additional  classrooms,  labs,  and  instructors  must  be  allocated. 
Larger  penalties  must  be  incurred  at  different  goal -intervals  as  the  various 
components  are  utilized  completely.  This  is  illustrated  graphically  by 
figure  2.  Certain  aspects  of  this  formulation  will  be  utilized  during  the 
development  of  the  specialized  algorithm  for  the  goal-interval  problem. 


3.  Development  of  Algorithm 

The  algorithm  developed  in  this  section  is  a  primal  simplex 
algorithm  which  will  exploit  the  specialized  structure  of  the  linear 
goal-interval  programming  problems.  Various  portions  of  the  algorithm 
can  be  recognized  as  modifications  and  extensions  of  other  techniques.  It 
is  due  to  the  ability  to  integrate  these  algorithms  into  the  development 
that  enables  computational  simplifications  that  will  provide  the  solution 
procedure.  Most  of  the  algorithms  and  techniques  employed  have  been 
mentioned  in  the  preceding  two  sections.  However,  there  is  one  additional 
algorithm  to  be  utilized  here.  It  was  constructed  by  Armstrong  and  Hultz  [3] 
to  study  restricted  case  of  the  Least  Absolute  Value  Estimator  (LAVE's) 


problem.  In  this  paper,  a  method  is  provided  for  a  “multiple  pivoting" 
technique  upon  which  some  extrinsic  characteristics  of  the  algorithm  of 
this  paper  will  rely. 

The  evolution  of  the  algorithm  will  begin  with  the  following  goal 
programming  problem  statement: 


subject  to 


m  p(k)  q(k)  _  . 

M.NIM1ZE  *  ■  111  *  l  W 
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+ 
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»p(k) 

(3.3) 

igk,i-i 

"  sM,k  =  *' 

2,  ... ,m;  i=l,2,... 

>q(k) 

(3.4) 

°i5ko 

-  9ko 

9ko’  k  =  2> 

•  •  • »  m. 

(3-5) 

It  should  be  noted  that  (3.1)-(3.5)  is  a  completely  general 
representation  of  the  formulations  of  the  problems  and  models  discussed  in 
section  2.  Any  linear  inequality  constraint  can  be  written  as  a  linear 
interval  constraint  of  the  form: 

bi  ~£1hijZj  ~  bi 

by  letting  b*  or  bT  (whichever  is  appropriate)  be  U  >  0  (or  possibly-U), 
where  U  is  a  non -Archimedean  transcendental  for  the  unbounded  case.  It  was 
stated  in  the  previous  section  that  these  types  of  constraints  were  due  to 
physical  limitations  on  the  system,  and  that  they  must  be  satisfied  before 
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the  attempt  is  made  to  resolve  the  problem  of  the  conflicting  goals.  Al¬ 
though  this  type  of  constraint  has  not  explicitly  incorporated  into  the 
system  (3.2)-(3.5),  the  situation  can  be  handle  by  assigning  an  arbitrarily 
large  positive  weight,  i.e.,  the  non-Archimedean  transcendental  M  >  0,  to 
the  deviations  associated  with  the  constraints  of  the  form  of  (3.6).  For 
simplicity,  it  is  assumed  that  the  system  associated  with  (3.6)  is  consistent. 

To  demonstrate  the  equivalence  between  the  goal  programming 
problem  given  by  (3.2)-(3.5)  and  the  goal  programming  problems  stated  in 
the  preceding  section  four  lemmas  are  stated.  The  proofs  for  these 
lemmas  are  extensions  of  work  found  in  the  references,  and,  therefore,  we 
shall  merely  indicate  the  appropriate  reference  rather  than  duplicate  results. 

LEMMA  1: 

Let  ( x ,6 )  be  an  optimal  solution  to  (3. l)-( 3. 5) .  It  then  follows  that 
P(k)^+  q(k) 

(  e  6  )  (  e  $k1)  =  0;  k  =  1,  2,  ...,  m. 

1  =  1  kl  1-1  kl 


PROOF: 

See  Chames,  Cooper  and  Ferguson  [11].  Note  the  assumption  that 


LEMMA  2: 


Let  ( x ,6 )  be  an  optimal  solution  to  (3.2)-(3.5).  It  then  follows  that 


M- 

6 


ki  "  9ki  '  gk,i-l  implieS  =  °  f0r  *  >  and  «ki  <  gk,i-l  -  9k, i  in,Plies 


\l  =  0  for  *  >  1- 


1*1  H  i  iraaMiMMnn  ilm-Jitr  i  iVT  m  ■  n  i  if  ■  ■  n  ■  «J 


PROOF 


See  Charnes  and  Lemke  [14].  Note  the  assumption  when  i  >  i . 

LEMMA  3: 

Let  (U)  be  an  optimal  solution  to  (3.1)-(3.5).  It  then  follows  that 

*k1  >  0  1mP'1es  Xkt  ■  4  -  4  .,(*«<!. 

and  6'ki  >  0  implies  2^  -  g^.j  -  g^for  t  <  i, 

PROOF: 

See  Charnes  and  Lemke  [14]. 

LEMMA  4: 

Let  r  be  the  rank  of  the  matrix  A  =  (a. .).  If  r  =  m  then  the  optimal 

ij 

objective  value  for  (3.1)-(3.5)  is  z  =  0. 

PROOF: 

See  Charnes,  Cooper,  Klingman  and  Niehaus  [12].  In  addition,  they  show  that 
if  r+l=m,  then  the  optimal  solution  can  be  obtained  explicitly. 

At  any  stage,  the  algorithm  will  work  with  a  submatrix  of  A 
which  has  full  row  rank  and  forms  a  basis  for  the  vector  space  spanned  by 

the  rows  of  A.  The  current  solution  will  be  such  that  a  basic  constraint 

holds  as  an  equality  at  one  of  the  extreme  points  for  the  goal  functional; 

that  is,  it  will  pass  through  g^  or  g^  for  some  i  whenever  k  is  basic. 

Also,  the  conditions  of  lemmas  1,  2,  and  3  are  satisfied  by  the  current 
solution.  The  fundamental  strategy  of  the  algorithm  is  as  follows: 

(i)  Without  varying  r-1  of  the  basic  constraints  (where  r  is  the 
rank  of  A),  determine  a  direction  of  change  for  the  remaining  basic  constraint 
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so  as  to  decrease  the  objective  function's  value.  If  it  is  not  possible  to 
obtain  a  decrease  for  any  of  the  basic  constraints,  then  the  current  solution 
is  optimal. 

(ii)  Assuming  that  the  value  of  the  objective  function  can  be 
decreased,  then  move  in  the  prescribed  direction  until  the  rate  of 
change  becomes  nonnegative.  Since  the  rate  of  change,  itself,  is  a 
piecewise  linear  convex  function,  then  it  will  have  changes  of  slope 
occurring  whenever  movement  along  the  constraint  intersects  an  extreme 
point  for  the  goal  functionals.  Therefore,  a  new  basis  will  be  formed  by 
the  previous  (r-1)  basic  constraints  which  remained  stationary,  and  the 
constraint  whose  value  caused  the  rate  of  change  of  the  objective  function  to 
become  nonnegative. 

Using  matrix  notation,  (3.1)-(3.5)  can  be  rewritten  in  the  following 


form: 


subject  to 


MINIMIZE  z  =  w+«+  +  W "6", 

(3.7) 

«  -  -  +  +  + 

(3.8) 

Ax  +  I  6  -  I  6  +IS  =  g 

0 

n  +  + 

0  <  6  <  Ag 

(3.9) 

0  <  6  <  Ag~ 

(3.10) 

0  1  6  1  A9g 

(3.11) 

where  all  vectors  and  matrices  are  defined  in  accordance  with  the  original 
summation  representation. 

The  matrix  A  may  be  partitioned  into  A  =  (^)  where  F  is  a  matrix 
with  full  row  rank.  The  algorithm  begins  by  choosing  such  a  matrix 
F  (see  Ben-Israel  and  Charnes  [6]  for  a  method  to  accomplish  this). 


The  complete  problem  is  now  rewritten  by  partitioning  it  into 


the  various  components  associated  with  the  matrices  F  and  R.  The  problem 
has  the  following  representation: 

MINIMIZE  2  =  w+/F  +  wRSR  +  w~6p  +  wr<5r 


subject  to 


«  ♦  IfS|: 


+  + 

Ipfip 


Rf#'  “  +  Vl  - 

„  +  + 

0  <  6  <_  Ag 


*F6F  «F6 


its*  +  Is  =  g 
R  R  R  R 


+ 

Ro 


0  <  S  <  fig 


0  <  5  <A  gQ 


where  a  =  Fx 

It  is  assumed  that,  at  any  stage  (iteration),  the  algorithm  has  a 

current  solution  which  satisfies  the  following  properties. 

(i)  The  conditions  of  lemmas  1,  2,  and  3  are  satisfied. 

(ii)  All  components  of  6p,  6p  and  6p  are  either  at  their  upper 
or  lower  bound. 

(iii)  For  each  nonbasic  goal  constraint,  exactly  one  representative 

+ 

component  of  6R,  6R  and  5R  is  between  its  upper  and  lower  limit. 

Property  (iii)  eliminates  the  difficulties  encountered  when  a 
degenerate  solution  exists.  These  difficulties  may  be  circumvented  by  a 
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perturbation  technique  similar  to  that  described  by  Charnes  [71.  To 
facilitate  the  algorithmic  description,  the  possibility  of  degeneracy  is 
ignored  here. 

An  initial  basic  solution  will  clearly  satisfy  (i)  and  (ii), 
and,  by  assumption,  it  will  also  satisfy  (iii).  The  following  develops  a 
procedure  for  moving  to  another  basic  solution  which  satisfies  these 
properties  and  decreases  the  value  of  z. 

At  this  point,  additional  notation  that  will  be  used  throughtout 
the  remainder  of  this  papaer  is  presented.  Let  be  the  component 

from  6p,6p,  and  6R  currently  not  at  its  upper  or  lower  bound. 

^  when  6^.  =  0  for  all  i ,  6ko  =  0 , 

6ki  ~  Agki  for  1  <  1  and 
0  <  6k«.  <  V  '  9k,Jl-l  =  A9kS, 

6ko  when  6ki  =  6ki  =  0  for  a11  1 

and  0  <  6ko  <  9ko  ■  9ko  =  A9ko 

when  6^  =  0  for  all  i, 

5ko  ~  Agko 

6m  =  for  1  <  * 
and  0  <  <  9k.e-l  -  9k,t  *  4% 
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A  notation  for  the  indexing  of  various  sets  of  also  required. 

IP  =  {k|A^x  >  g*p  and  Ak  is  a  row  of  R} 

IN  =  { k |  Akx  <  gkQ  and  Ak  is  a  row  of  R} 

IM  =  { k| g“Q  <  Akx  <  gkQ  and  Ak  is  a  row  of  R> 


and  IB  =  { k| A^  is  a  row  of  F  or  Ak  is  a  basic  row}. 

The  algorithm  will  determine  if  a  variation  in  a  particular 
component  of  a,  say  ap,  will  decrease  the  value  of  the  objective  function. 
All  other  a.,  for  i^p,  will  remain  at  their  current  value.  This  perturbed 

A  A 

vector  is  denoted  by  a.  Its  components  will  be  cu  =  ou  for  i^p,  and 

a  =  a  +  e.  If  e  is  chosen  sufficiently  small  so  that  none  of  the 
P  P 

elements  will  move  to  their  upper  or  lower  bound,  the  rate  of  change 
of  the  objective  function  will  be 


+  Sp  wk,s(k)AkF*p  "  iSn  wk,s(k)AkF*p)si9n(e) 


w* 


IB  ( p ) 


Vzp(e)  =  (wIB(p) 

WIB(p),£+l  if  ap  =  9k£  E  9IB(p),Jl  3nd  e  *  0 
+  . 

wIB(p) ,Jl  if  ap  =  gk&  and  e  <  0 

and  e  >  0 


wIB(p) ,i  if  ap  =  g 


k£ 


wIB(p) ,£+l  if  ap  =  gkJl 


and  e  <  U 


where  we  define  WjB{p)>Q  =  wJB(p)>Q  =  0 

M  # 

and  F._  is  the  p-th  column  of  F  and 
P 


sign(e)  = 


e  >  0 
e  <  0 
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If  the  gradient  Vzp(e)  is  negative,  determine  how  ap  can  be 

A 

varied  before  the  rate  of  change  becomes  nonnegative.  As  ap  varies,  the 
gradient  Vzp(e)  will  change  whenever  any  reaches  an  upper  or 

lower  bound,  and,  therefore,  it  is  necessary  to  recalculate  Vzp(c)  at  each 
of  these  new  extreme  points.  This  procedure  of  determining  the  new 
gradients  will  be  further  elaborated  upon  later  in  this  section.  At  this 
point,  the  main  aim  is  establish  a  procedure  to  find  the  sequence  of  the 
extreme  points.  The  procedure  is  an  extension  of  the  standard  minimum 
ratio  test.  From  the  system  of  nonbasic  constraints, 

RF  a  +  IR6R  -  Ir6r  +  IR6R  =  gRQ,  it  is  recalled  that  there  will  be  only 
one  component,  for  each  constraint  which  is  not  at  its  upper 

or  lower  bound. 

The  first  change  in  Vzp(e)  will  occur  for  the  minimum  of  the 
following  ratios.  The  possible  ratios  are  broken  down  into  various 
cases  and  a  brief  explanation  for  each  case  is  provided. 


(1) 


(2) 


0*{A) 


Minimum* 
k  ( 


Ag* 


Ki.s(k) 


6k,s(k) 


l\F* 


if 

;  sign(e)Aj(F.p  >  0  and 

k  e  IP  or 
£ 

sign(t)AkF"p  <  0  and 
k  e  INU  IM 


increases  to  its  upper  bound 

*  (6k. s(k)  # 

4>  (A)  =  M1nimum< - 3 —  ;  sign(e)A.F. 

k  l|AkF*pl  k  * 


<  0  and 


k  c  IP  or 

sign(E)AkF?p  >  0  and 
k  e  INU  IM 


6*  decreases  to  zero 


■9!b(p).»+1  -  9fB(p).lt<  •  £  ’  0  "“  *  *  * 

e  <  0  and  *  *  - 

I»Tb(p),1  -  «!b(p).»-iI  £<0  "d*** 

£  >  0  and  *  *  - 

where  we  define  gJB(p)>.j  •  gIB(p)t0  *"d  g^^.j  -  g;B(p)>0 

This  Is  the  possible  change  in  oip  before  £  reaches  an  upper  bound 

or  fi|B(p),£-l  reactes  a  ^ower  bound.  Therefore,  the  first  extreme  point 
along  the  path  will  be  associated  with  the  minimum  ratio: 

0  «  MIN  (0*(A)b4*(A).O*) 

The  algorithm  now  updates  the  gradient  for  this  new  extreme  point.  The 

description  of  this  process  is  also  broken  down  Into  the  cases  which 

correspond  to  those  of  the  minimum  ratio.  The  following  Indicates  the  term 

to  be  added  to  Vz  (e). 

P 


(1) 

K,s(kH--:.s(k)>i*/pi 

u  e  IP 

(2) 

("I,s(k)  •  "!,,(k)-i>lA.F-pl 

(3) 

W!b(p),£+1  '  W!b(p),£ 

and 


“ie(p),i  "  *!«(P).i-i 

where  we  define  wfc0  ■  wk'0  -  w*0  ■  0,  w*_.j  -  wfcl  -  wjjj 
and  "k,-l  '  wkl 


where  u  corresponds  to  the  row  yielding  the  minimum  ratio. 


When  the  new  gradient  has  been  calculated,  then  we  shall  either 
proceed  to  the  next  extreme  point  on  the  path  it  is  has  remained  negative, 
or  we  shall  perform  the  standard  simplex  pivot  operation  if  it  has  become 


nonnegative.  If  the  pivot  is  performed,  then  IB(p)  will  leave  the  basis 
and  be  replaced  by  the  equation  in  R  associated  with  the  minimum  ratio. 

If  the  new  rate  of  change  is  still  negative,  the  algorithm 
updates  the  particular  ratio  which  provided  the  minimum.  All  other  ratios 
will  remain  the  same.  To  provide  simplicity  in  the  exposition  and  also  to 
illuminate  the  process,  we  shall  assume  that  the  minimum  ratio  was 
associated  with  case  (1)  and  that  row  r  yielded  this  ratio.  Thus,  the 
distance  to  the  next  extreme  point  (or  goal-level)  for  the  equation 
associated  with  this  minimum  ratio  can  be  found  by  adding  the  width  of  the 
next  goal -segment  to  the  numerator  of  the  ratio.  This  new  ratio  is  denoted 
as  follows: 


0  (A+Ax)  = 


K«.  - 8 


T+  \  .  + 

rV  9r1q(r) 


i\PF' 


sign(e)*AkpF#  >  0 


The  new  ratio  0+(A+Aj)  replaces  the  ratio  in  the  r-th  row,  and  the  new 
minimum  ratio  is  determined.  The  rate  of  change  is  then  updated.  This 
process  is  repeated  until  the  rate  of  change  becomes  nonnegative,  and  at 
that  point  the  pivot  procedure  is  implemented,  if  the  minimum  ratio  is  not 
given  by  0^. 


An  additional  aspect  of  the  algorithm  is  that  some  bookkeeping 
is  required  to  accomodate  the  movement  of  6k*  for  cases  (2)  and  (4), 
respectively.  If  the  multiple  pivoting  technique  proceeds  along  an 
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extreme  point  path  starting  at  <sj^,  for  example,  then  the  variables  associa¬ 
ted  with  the  decreasing  sequence  of  segments  would  be  5k(£-2) * . 

and,  6^.  If  further  decreases  in  the  objective  function  are  possible, 
then  pass  into  the  goal -Interval .  Thus,  as  the  algorithm  through  the  various 
types  of  intervals  associated  with  the  index  sets  IP,  IN,  and  IM,  this  pro¬ 
cess  must  be  monitored. 

Application  of  this  pivoting  procedure  will  result  in  the  optimal  sol¬ 
ution  being  obtained  in  a  finite  number  of  iterations  (degeneracy  being  re¬ 
solved  in  an  appropriate  manner).  The  algorithm  increases,  or  decreases, 
the  value  of  a  "basic"  goal  constraint  until  any  further  unilateral  change 
would  increase  the  value  of  the  goal  functional.  The  end  result  of  this 
process  is  to,  at  times,  combine  several  standard  iterations  into  a  single 
one.  Thus,  basis  updates  and  the  recalculation  of  ratios  are  reduced. 

4.  Conclusions 

Goal  programming  has  long  been  recognized  as  one  of  the  most  popular 
modeling  techinques  in  operations  research.  Despite  the  many  applications 
of  goal  programming,  little  research  seems  to  have  been  directed  towards 
special  algorithms  to  solve  the  resulting  problems.  One  recent  effort  in 
this  area  is  given  by  Charnes,  Cooper,  Klingman  and  Niehaus  [12].  In  this 
paper,  we  have  developed  a  specialized  algorithm  to  solve  goal-interval 
programming  problems.  The  algorithm  is  able  to  exploit  the  unique  structure 
of  problems  formulated  in  this  manner.  It  possesses  a  capability  to  combine 
several  standard  simplex  pivots  into  a  single  pivot.  Computational  studies 
with  various  computer  code  implementations  of  the  algorithm  will  be  the 
subject  of  future  research. 


j 
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While  the  main  thrust  of  this  paper  has  been  devoted  to  algorithmic 
considerations,  the  importance  of  the  model  itself  should  not  be  overlooked. 
Many  goal  programming  problems  found  in  the  literature  can  be  expressed  in 
the  framework  of  (3. l)-(3.5).  Refinement  of  the  goal  intervals  will  also 
allows  us  to  obtain  an  c-optimal  solution  to  a  separable  nonlinear  goal 
functional.  An  immediate  consequence  of  this  is  to  allow  consideration 
of  the  L  norm  1  <  p  <  «>. 

p  _  r  _ 
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